Spatial entanglement in optical parametric oscillators with photonic crystals 
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The effects of an intracavity photonic crystal in a multimode optical parametric oscillator are 
studied, with a special focus on quantum fluctuations. The capability to either stimulate or inhibit 
the spatial instability, lowering or increasing the parametric threshold, allows to control the inten- 
sity fluctuations and correlations. A significative quantum noise reduction and an increase of the 
range of squeezed quadratures are found above threshold where spatial Einstein-Podolsky-Rosen 
entanglement and inseparability are found. 
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Photonic crystals (PC) are dielectric media with peri- 
odic modulation of the refractive index which can lead 
to gaps in the allowed frequencies of electromagnetic 
waves [1] . Seminal works predicted the possibility to con- 
trol spontaneous emission by PC [2 when suppressing a 
radiative transition with frequency within the photonic 
band-gap. After two decades of intense experimental ac- 
tivity, these engineered media provide an unprecedented 
control of light confinement, guiding, and propagation 
[U |3] . Recently, inhibition of spontaneous emission |4] 
has been experimentally shown, improving the extraction 
efficiency of light emitting devices [5] and redistributing 
where needed the corresponding energy [6]. The use of 
PC for environment (dissipation) engineering is also the 
basis of intense research activity about non-Markovian 
evolution of quantum states j7] . An unexplored issue we 
raise here is if a PC can be used to control or even sup- 
press light fluctuations to improve quantum correlations 
in the continuous variable regime [8 in driven devices 
emitting entangled beams. Quantum correlated bright 
beams can be generated by optical parametric oscilla- 
tors (OPOs). These are common devices giving squeez- 
ing [9] and entanglement ^ between modes with dif- 
ferent polarization or frequency. In 1992 similar quan- 
tum effects were also predicted between different spatial 
modes extending the study of quantum effects to 

multimode OPO [HIIS], Kerr media [II1[14], and second 
harmonic generation [15] . Different applications based on 
spatially multimode operation have been recently demon- 
strated in optical switching [16], quantum imaging [17 , 
metrology [18], and quantum information [19j. Consid- 
ering both theoretical activity and experimental achieve- 
ments, the subject of spatial entanglement is becoming a 
mature research topic. 

In this Letter we explore how to control and improve 
multimode squeezing and entanglement by means of an 
intracavity photonic crystal (PC) in a nonlinear device. 
The prototype system considered is a multimode degen- 
erate OPO, well studied both below and above thresh- 
old [T2I [T^, when a PC is introduced. Within the 



proposed photonic-crystal optical parametric oscillator 
(PCOPO), the spatial modulation is not changing the 
environment spectrum of fluctuations [7^ but is instead 
modifying the Hamiltonian of the intracavity process. 
The PC modulation takes place in the transverse plane 
-the device is longitudinally monomode- and is mod- 
eled by a spatial profile of the otherwise homogeneous 
refractive index. Intracavity PC have been predicted to 
allow for inhibition of the phenomenon of pattern for- 
mation in Kerr and singly resonant OPOs [23l [24] by an 
increase of the pump energy needed to cross the insta- 
bility threshold. The phenomenon arises when unstable 
spatial modes are in the band-gap because their emis- 
sion is prevented, as recently observed in two indepen- 
dent experimental set-ups [25]. Here we will show that 
in PCOPO the possibility to tune the spatial instability 
- matching-up to the parametric threshold- is actually 
wider. Indeed, the PC allows to either increase or re- 
duce the threshold energy due to the mixing of differ- 
ent frequency waves. Related to the instability shift is 
the possibility to control the fluorescence intensity below 
threshold and the spatial distribution of quantum fluc- 
tuations. After characterizing quantum correlations in 
different regimes, two major effects are identified above 
threshold. First of all, the PC provides a locking mech- 
anism to freeze the known diffusive drift motion of the 
pattern \W reducing by orders of magnitude the fluctu- 
ations in the associated quadrature. A larger range of 
quadratures is then found to be squeezed, decreasing the 
sensitivity to the choice of the local oscillator phase in 
squeezing measurements. Moreover, we find that spatial 
modes can be entangled above threshold when introduc- 
ing the PC in the OPO, either considering state insepa- 
rability [21] or the Einstein-Podolsky-Rosen criterion of 
Ref. [^. 

The master equation for multimode type I degenerate 
OPO was described in detail by Gatti et al. [12 . A de- 
scription valid both below and above threshold (for pump 
|ao| < 2) can be obtained through a mapping into the 
Q-representation as discussed in Ref. [13]. This leads to 
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nonlinear Langevin equations for the spatially dependent 
pump ao and signal ai fields 

dtao{x,t) = -[(l + iAo(x))-iV^]ao(x,t) + 

dtai{x,t) = -[(l + zAi(x))-2zV^]ai(x,t) + 

ao{x,t)al{x,t) (1) 

with E input field, diffraction, and additive and 
multiplicative phase-sensitive white noises [28 . The PC 
refractive index modulation is modeled by introducing 
spatial dependent detunings Ao{x) and Ai(x), which can 
have different amplitudes and, in the simplest case, have 
the same periodicity with wave-number kpc. The main 
mechanism we aim to explore is the effect of the band- 
gap on the spatially multimode down- conversion process. 
As known, modulation instability in OPO with negative 
signal detuning appears at wave-number kc = y^— Ai/2 
^26|. Therefore, the most interesting configuration is for 
a (sinusoidal) modulation with kpc = 2kc^ since in this 
case the signal pattern would be in the photonic band- 

gap m- 

We start considering the effects on the quantum fiuctu- 
ations in the PC OPO below threshold. On average the 
signal field vanishes everywhere = but, due 

to the nonlinearity of the medium, it is not in a coher- 
ent vacuum state, either with or without the PC. When 
spatial modes are considered in the far field, the ampli- 
tude of intensity fiuctuations is maximal at the critical 
wavenumber [12] (Fig. [ij)). As the least damped spa- 
tial mode falls within the PC band-gap, inhibition of the 
off-axis emission is expected [11 [4]. This effect can be ap- 
preciated comparing the larger intensity fiuctuations for 
the OPO (Fig. [T]3, continuous line) with respect to the 
PCOPO with modulation of the signal detuning Ai(x) 
(Fig.[l]3, dashed line). On the other hand, it is rather sur- 
prising to find that, for a fixed input energy, fiuctuations 
in the PCOPO can also be increased by the modulation 
of the refractive index introduced by the PC (Fig. [l] dot- 
dashed lines). In other words, when the signal detuning 
Ai is modulated, the intensity fiuctuations are lowered 
[4], while when (also) the pump detuning Aq is modu- 
lated, such fiuctuations increase, in spite of being in the 
PC band-gap. 

The possibility to either lower or increase the fiuctua- 
tions of the most intense spatial modes is due to the pres- 
ence of wave-mixing between different frequencies in the 
parametric oscillator. Indeed, the fiuctuations strength 
is inherently related to the proximity to the instability 
threshold in a nonlinear system driven out of equilibrium 
and approaching this point. Large fiuctuations are a clear 
signature of an instability and this allows a consistent 
description of our results. In Fig. [l^, analytical average 
intensities in the linear approximation are shown for 
the PCOPO below threshold, in different configurations. 



When the threshold is raised (pattern inhibition due to 
the PC [23l |24]), the fiuctuations strength is reduced, 
while -for the same input energy- if the threshold is low- 
ered by introducing the PC then fiuctuations increase. 
What we actually find is that, in presence of nonlinear 
mixing between different fields, the way in which the PC 
changes the instability of the signal is not trivial. We 
remind that when the stripe pattern at the critical kc ap- 
pears in the signal field, the pump also develop a modula- 
tion at 2kc' Then for the detuning modulation k^c = 2kc 
the unstable signal wavenumber is in the band-gap but 
the detuning profile somehow 'stimulates' the nonlinear ly 
emerging pattern in the pump field. Therefore this PC 
modulation provides two competing mechanisms, inhibit- 
ing the signal spatial instability, as in Refs. [23l[24], but 
also imprinting in the pump the nonlinear structure fa- 
voring the instability process. The vertical asymptotes 
in Fig. [l^, where the linear approximation for analytical 
calculations breaks down, show the different threshold 
values for E depending on the PCOPO configuration. 
Notably, if the PC modulates the pump detuning, the 
parametric threshold can be crossed for values even lower 
than in the case of perfectly resonant OPO. 
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FIG. 1: (Color online)(a) Steady intensity {A\{kc)Ai(kc)) in- 
creasing with E (below threshold). Here the symbols are re- 
sults of numerical simulations of Eqs. ([T]) at different pump 
values, while lines represent analytical results within a linear 
approximation [27]. (b) {A\{k)Ai{k)) from numerical sim- 
ulations for E = 0.9. Different curves correspond to the 
OPO without PC with Aq = and Ai = -1 (black solid 
line); Ai = — 1 + 0.5 sin(/cpca:) and Ao = (green dashed 
line); Ao = 0.5 sin(A:pcx) and Ai = — 1 (red dot-dashed line); 
Ao = 0.5 sin(A;pca:) and Ai = — 1 + 0.5 sin(/cpca:) (orange 3 
dots-dashed line). 

Apart from the strength of spatial fiuctuations, an im- 
portant aspect is the quantumness of the correlations. 
Non-classical effects in multimode OPO are known to 
exist between opposite far field modes +/c and —k due 
to emission of photons pairs in the parametric down- 
conversion process [TTJ [12]. In particular, two- modes 
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squeezing is studied considering the generic joint quadra- 
ture 

^0^{k, -k) = {Ai{k) + XAi{-k)e'^)e'^ + h.c. (2) 

with (/) relative phase between the superposed spatial 
modes. Here we take A = 1. Squeezing achieved be- 




FIG. 2: (Color online) Variance of Sg,^ (from Eq. p| with 
A = 1) as a function of the quadrature angle and for the su- 
perposition angle 4> giving the largest squeezing for each OPO 
and PCOPO configuration. Pump field 5% below threshold 
(a) and 2% above threshold (b). The horizontal dashed line 
is the shot noise and other lines as in Fig.[l] The smaU arrow 
shows the deviation between different numerical runs in the 
OPO, due to phase diffusion. 

low threshold increases with the pump intensity being 
maximum at the parametric threshold pTHT3] . Due to 
the discussed PC effect on the parametric threshold, the 
squeezing attained in OPO and PCOPO will be com- 
pared at the same distance from the respective thresh- 
olds. We find then that squeezing achieves similar values 
in the OPO and in the PCOPO modulating pump and/or 
signal detuning, as shown in Fig. [2^, being the major dif- 
ference in the dependence on the angles and (j). 

A different scenario is found above threshold, consid- 
ering squeezing between intense modes. In Fig. ^ the 
variance of Sg*^ (Eq. |2| for the OPO (upper line) is com- 
pared with the PCOPO (three lower lines). Even if the 
attained squeezing (minimum value of the plotted vari- 
ance) is similar in all cases, there are important differ- 
ences in the noise present in the unsqueezed quadrature 
(maximum value). Far from being in a minimum uncer- 
tainty state, the OPO displays an extremely large noise in 
the unsqueezed quadrature (black line in Fig. ^p) due to 
the well-known phase diffusion between down-converted 
modes and to excess noise in their relative phase [20] . In 
spatially multimode devices, the visible effect is a diffu- 
sive motion of the excited pattern and has been related 
to translational symmetry break and noise excitation of 
the corresponding neutral Goldstone mode [14 . On the 
other hand, due to the refractive index modulation, the 
PCOPO does not exhibits translational symmetry and 
the formed pattern is locked to the position of the PC. 
This leads to a strong reduction of noise in the PCOPO: 
in Fig. |2]d there are two orders of magnitude between the 
variances of the unsqueezed quadrature in the PCOPO 
for modulated Aq and the OPO, as highlighted by the 



horizontal gray stripe. An important consequence in 
view of applications is that the reduction of fluctuations 
in unsqueezed quadratures leads to a significant increase 
of the range of quadratures with sub-shot-noise fluctua- 
tions, as highlighted by the two vertical stripes in Fig.[2]3. 
A PCOPO will indeed be more robust to changes in the 
choice of the local oscillator phase 6 as more quadratures 
are actually squeezed. 




FIG. 3: Variances sum X (a and b) and product S (c and 
d), as defined in Eqs. Q and ^ for the output fields [30] 
of OPO (a and c) and of PCOPO (b and d) when varying 
the superposition and interference angles , 0o • In (b) and 
(d) the PC modulates the pump, Ao = 0.5sm{kpcx) (like 
for the red dot-dashed line in Fig. [l] and [2|) and all results 
are for pump 2% above threshold. The white dashed lines 
limit wide regions regions for which EPR entanglement and 
inseparability are respectively predicted. 

Multimode OPOs allow to generate not only squeezed 
but also spatially entangled states [El A series of 
key experiments recently demonstrated spatial entangle- 
ment between light beams (continuous variable regime) 
in different optical nonlinear devices [17 . Here we show 
how entanglement in parametric oscillators is changed 
by the presence of a photonic crystal, considering two 
well-known criteria [20] [21] . One distinguishes states ex- 
hibiting Einstein- Podolsky- Rosen (EPR) paradox [31 for 
conditional variances such that 

£ = A^i;0o,0oA^i;^o+-,0o+7, < 1, (3) 

for some choice of superposition and interference angles 
^0,00 l2Qi- Notice that here the parameter A minimizes 
each variance of the joint quadrature as also reviewed 
in Ref. [l^. A second measure we consider is the insep- 
arability condition 

I = A^S,,,^, + A2S,,+5,^,+^ < 2{a' + ^). (4) 

with = {aAi{k) +a"^Ai(-fc)e*^)e*^ + /i.c. and pos- 
itive parameter a [21 . Below threshold, the presence 
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of the PC changes significantly the intermode correla- 
tions of this systems leading to new not vanishing terms 
with respect to the case of the OPO (for instance in the 
PCOPO {Al{k)) 7^ 0), but as for squeezing, similar re- 
sults are found for the attained entanglement. The most 
significant differences are found, again, above threshold 
where the presence of the PC enhances quantum effects 
leading to a spatially entangled state. In this regime, 
the mentioned phase diffusion [Mj |20l leads to spikes at 
low frequency in noise spectra, preventing entanglement 
in the OPO, as we show in Fig. [3^ and c. On the other 
hand, for a non translational invariant system such as the 
PCOPO, we find significant regions in which both EPR 
paradox (|3| and state inseparability Q are predicted. 
We stress that in Fig.jSj OPO and PCOPO are compared 
at the same distance form the instability threshold and 
that, even if the best performance is obtained when pump 
detuning is modulated (Fig. [sJd and d) we find that all 
configurations of the PCOPO show some entanglement, 
degraded when removing the PC. 

Summarizing, we have analyzed the effect of a PC in 
a multimode OPO whose detuning suffers a transverse 
spatial modulation, in the case in which the spatially 
unstable mode appearing at threshold falls within the 
band-gap. Due to the presence of nonlinear mixing be- 
tween waves at 2uj and uj there are competing phenomena 
leading to either inhibition and enhancement of intensity 
quantum fluctuations for a fixed pump below threshold. 
The possibility to control with the PC the intensity of 
quantum fluctuations is related to the raise (pattern in- 
hibition) and lowering of the parametric threshold. We 
find that, below threshold, the attained spatial squeezing 
as well as entanglement in the signal field are preserved at 
a fixed distance form the parametric threshold. Notably, 
above threshold, the break of the translational invariance 
due to the PC provides a strong mechanism to reduce (up 
to two orders of magnitude) the quantum fluctuations 
associated to pattern diffusion leading to squeezing over 
a significantly larger range of quadrature angles. This 
would reduce the sensitivity of the measurable squeezing 
with the choice of the phase of the local oscillator. More- 
over, the strong spatial locking due to the presence of the 
PC in the OPO allows to generate inseparable as well as 
EPR entangled spatial beams. This analysis, restricted 
to one transverse dimension, is expected to give similar 
results when extended to 2D for PC modulated only in 
one direction, as also confirmed by our analytical results 
below threshold [27 . More complex is the case of a PC 
with transverse hexagonal or square geometries where the 
same pattern selection process is an open question. 
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